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Abstract 

ON 

We study the neutrino mass matrix in supersymmetric models in which the 



quark and charged lepton mass hierarchies and also the suppression of baryon 
or lepton number violating couplings are all explained by horizontal U(l)x 



^ ■ symmetry. It is found that the neutrino masses and mixing angles suggested 



by recent atmospheric and solar neutrino experiments arise naturally in this 
framework which fits in best with gauge-mediated supersymmetry breaking 
with large tan /3. This framework highly favors the small angle MSW oscil- 
lation of solar neutrinos, and determine the order of magnitudes of all the 

neutrino mixing angles and mass hierarchies. 
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Typeset using REVIgX 



The fermion mass problem consists of understanding the flavor mixing structure among 
quarks or leptons as well as the hierarchy of their mass eigenvalues. It has been suggested 
that these hierarchical structures can be explained by a horizontal U(l)x symmetry whose 
spontaneous breaking is described by A ~ Cabbibo angle |1|-§J. Recent experimental data 
on atmospheric and solar neutrinos suggest non- vanishing neutrino masses and mixing . 
If spontaneously broken U(l)x is the origin of the quark and charged lepton mass spectrum, 
it is expected to have implications for the neutrino masses also. It has been noted that when 
implemented in supersymmetric (SUSY) models, U(l)x can explain not only the quark and 
lepton mass spectrum, but also the smallness of dangerous baryon/lepton number (B/L) 
violating interactions Q. This framework is interesting since renormalizable L- violating 
couplings are small enough to satisfy the current experimental bounds, but still nonvanishing 
and thus can generate neutrino masses. In this paper, we wish to examine the possibility 
that the neutrino masses and mixing angles suggested by recent atmospheric and solar 
neutrino experiments arise naturally in the framework of SUSY models in which the quark 
and charged lepton mass hierarchies and also the suppression of B/L- violating couplings 
are all explained by horizontal U(l)x symmetry. Combining the neutrino oscillation data 
with the informations from the quark and charged lepton sector and also the constraints on 
B/L- violating couplings, we find the U(l)x charge assignments producing all the fermion 
masses and mixing angles correctly. This framework fits in best with gauge-mediated SUSY 
breaking models with large tan (3, favors the small angle MSW oscillation of solar neutrinos 
over the large angle just-so oscillation, and determines the order of magnitudes of all the 
neutrino mixing angles and mass eigenvalues. In this framework, 11x2/ vn^ m 4 x 10~ 2 is 
essentially due to the loop to tree mass ratio, while mi/m 2 ~ U^ 2 ~ A 4 is due to the U{l)x 
selection rule where (A = 1,2,3) denote the neutrino mass eigenvalues and UiA the 
mixing matrix. 

To proceed, let us briefly summarize the U(l)x selection rule estimating the size of 
couplings [T]-|3| . The Kahler potential and superpotential of the model are generically given 
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by 



K = Z u (\, A)$ 7 <T J + [X/j(A, A)$ / $ J + h.c] + 
W = ^/i/j(A)$ 7 $ J + ^F W (A)$ 7 $ J $^ + .., 



(1) 



where $ J denote light chiral superfields and the ellipses stand for the terms of higher order 
in & . The £/(l)x-breaking order parameter A corresponds to the VEV of a chiral superfield 
4> with the U(l)x charge X(<p) = —1: A = {4>)/M « 1/5 for the fundamental mass scale 
M which is presumed to be of order the Planck scale M P . The U(l)x selection rule states 
that the hierarchical structures among the coefficients are due to the insertion of A = ((f)) /M 
or of A = ((f)*) /M to make the corresponding operators to be [/(l)x-invariant. This leads 
to Z u « \\*'- x j\, X u « Al^+^l , flu « fx\ Xl+Xj - x W, Y IJK « \ x i +x -> +x Ke( Xl + xj + x K ) 
where xj = X( ( l> / ), i.e. the U(l)x charge of $ / , jl denotes the representative component 
of flu whose operator has the U(l)x charge X(p), and 6(x) = 1 when a; is a non-negative 
integer, while 9(x) = otherwise. The overall size of dimensionful flu depends upon the 
mechanism generating the corresponding bilinear terms and can differ from the fundamental 
mass scale M in general. 

After integrating out supersymmetry breaking fields while taking into account supergrav- 
ity effects, one can redefine the chiral superfields in the resulting effective theory to have a 
canonical Kahler metric: $ 7 — > Ru& J where R L j obeys RjjZ jk R* kl = 5 IL . The order of 
magnitude estimate of Zu above implies Ru \\ x i- x j\. Then for the redefined $ / with 
canonical kinetic term, the bilinear and trilinear couplings of the effective superpotential are 



including first the contribution from the bare superpotential W, second the effects of super- 
field redefinition $ 7 — > i?/j$ J , and finally the supergravity contribution from the Kahler 



given by 



jW/j ~ p>u + RikRjl&kl + m 3/2 Xu 



Yijk ~ Yijk + RilRjmRknYlmn, 



(2) 



potential which is proportional to the gravitino mass m 3 / 2 . 
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The most general SU(3) C x SU(2)l X [/(l)y-invariant superpotential of the MSSM su- 
perfields is given by 

+A% k U?DjD° k + Xy.^ir, + KikULjEi 

where (Hi,H 2 ), {Li, El), (Qi,U?,Df) denote the Higgs, lepton, and quark superfields, re- 
spectively. Among the possible non-renormalizable operators, we include only the d = 5 
see-saw operator which is presumed to be induced by the physics at the scale Ms- Here 
we have rotated away the possible bilinear term /ijLjif 2 through the unitary rotation of 
superfields, which does not alter the order of magnitude estimates of couplings. 

The U{l)x charges denoted by the small letters qi,Ui, e.t.c. for the superfields Qi,U?, 
e.t.c. are well constrained by the experimental data. The quark Yukawa couplings Y^j' 
are determined by the U(l)x charges of the operators H 2 QiUj and H\QiD c -. The large 
m t « (H 2 ) suggests first of all q 3 + u 3 + h 2 = 0, while tan (3 ~ \ x m t where x = q 3 + d 3 + hi 
and tan f3 = (H 2 ) / (Hi). Then the observed pattern of the CKM matrix and the quark mass 
eigenvalues (m u ) ~ m 4 (A 8 ,A 4 , 1) and (m d ) « mh(A 4 ,A 2 ,l) contain enough informations to 
reconstruct the Yukawa matrices Y u and Y d , leading to 

(?13, 923,^13,^23, d i3 ,d 23 ) 

= (3,2,5,2,1,0) or (-3,2,11,2,7,0), (4) 

where qij = qi — qj, e.t.c, and x < 3 for the second case. The charged lepton mass hierarchy 
(m e ) « m T (A 5 , A 2 , 1) suggests (ei 3 , e 23 ) = (5 - h 3 , 2 - 1 23 ) (Case I) or (9 - 1 13 , -2 - 1 23 ) (Case 
II) where = e, — ej, e.t.c. In order to have Y^i ~ ^ 2+x an d a ls° a non-singular Y2, 
we need more constraints: h 3 > l 23 or l\ 3 < l 23 — {x + 2) for Case I, —2 < Z 2 3 < 0, x < 2, 
Z13 > l 23 + 4 for l 23 > x — 2 for Case II. As we will see, the neutrino oscillation data imply 
^23 = and l\ 3 = 2, which is compatible only with Case I. We thus have 
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(ei3,e 23 ) = (5-Z 13 ,2-Z 23 ) = (3, 2). 



(5) 



The lepton flavor mixing is usually read off from the non-trivial neutrino mass matrix 
rn'ij in the basis of the diagonal charged lepton mass matrix. Denoting the unitary mixing 
matrix by U^, the mass eigenvalues are given by = Ui/Jj r jA m ij (A = 1,2,3). Recent 
Super-Kamiokande and other experiments on the atmospheric and solar neutrino oscillations 
HH suggest the following neutrino oscillation parameters: 



l 32 



Am-- ~~~ 



Am 2 31 « 2.2 x l(T 3 eV 2 



9 



~ 1 

atm ~ J- ? 



Am^ « 



5 x lCT 6 eV 2 , 6 sol « 3.7 x 1CT 2 



6.5 x lCr n eV 2 



7 sol 



0.52 



<- MSW 
just-so 



(6) 



where Am 2 



AS 



m 2 ^ — m^, 6> atm = 6* 32 or and # S oi = 9^1 or ^ for the effective mixing 



angle: sin 2 29% = A\U ai U* aj U* pi U Pj \. 

In our framework, L-violating couplings are suppressed by having Z, significantly bigger 
than hi and —h 2 . As will be shown explicitly later, the covariance under U(l)x suggests 
that the neutrino mass matrix in our framework takes the form: 



(m")« = m 3 X li3+lj3 A 



(7) 



where m 3 is the largest mass eigenvalue, Zj 3 = Zj — Z 3 , and all Ay are of order unity. This form 
of m 1 " leads first of all to Z7j 3 ~ X la . The large atmospheric u^-u r mixing unambiguously 
implies ~ 1, and thus 



z 2 = z,. 



(8) 



For Z 2 = Z 3 , the mass matrix (0) implies also Z7j 2 ~ A^ 3 . Combined with the unitarity, this 
determines the mixing matrix to take the form: 



U 



1 1 A /l3 A' 13 ^ 
A' 13 1 1 
A' 13 1 1 



(9) 



Given the above form of U, the MSW mixing angle 8 so i pa A 2 implies 

k = l 2 +2 = l 3 + 2, (10) 

while the just-so oscillation leads to h = l 2 = h- Recent Super-Kamiokande and CHOOZ 
data 0,0] indicates that rarely if ever oscillates into u e , which can be interpreted as 
excluding h = l 2 = h [H- 

We have seen that all Uia resulting from the mass matrix (0) are determined essentially 
by the U(l)x charges U. As it will become clear later, although all A*j are of order unity, the 
corresponding matrix is naturally approximately singular and thus gives a mass hierarchy 
m 3 ^> m 2 ^> mi. Then the oscillation data (||) implies 

m 3 ^5xl(T 2 eV, m 2 /m 3 pa 4 x 1(T 2 . (11) 

Let us now discuss how the neutrino mass matrix (|7]) with approximately singular Aij 
arises in SUSY models with U(l)x- Although not a unique possibility, an attractive scheme 
to suppress dangerous L- violating couplings in our framework is to have Zj significantly bigger 
than hi and —h 2 . In this scheme, one can easily arrange the physics at Mg, e.g. the U(l)x 
charges of the superheavy singlet neutrinos, to make the resulting see-saw coefficients 
in © suppressed by \ l ^+h+ 2h ^ ||. if M s is the string scale M str i ng pa 5 x 10 17 GeV or the 
unification scale Mqut ~2x 10 16 GeV, which is perhaps the most plausible possibility, this 
would result in 

m seesaw _ _ 1Q -4) x X k+l j+ 2h 2 e y ? ^ 

which is too small to be relevant for the atmospheric and solar neutrino masses for Zj signif- 
icantly bigger than —h 2 . In fact, the two representative models that we found in this paper 
have li + h 2 > 8 and thus a completely negligible see-saw contribution. 

Once the see-saw contribution is negligibly small, the atmospheric and solar neutrino 
masses arise from the renormalizable interactions in the superpotential (0) and also the 



following soft SUSY breaking terms [J1Q| pT^I : 



(13) 



where now all field variables denote the scalar components of the corresponding superfields. 
The L-violating Bi or m\. Hl (in the basis where [liLiHi in the superpotential are rotated 



away) results in the tree-level neutrino mass [|13l|14 



m 



tree 



(14) 



where M a denote the 577(2) x U(l) gaugino masses and the sneutrino VEV's are given by 

2M z {m 2 L . Hl cos (3 + Bi sin (3) 



m~ 



±Mf cos2/3 



(15) 



for the Z-boson mass Mz and m ; ~ denoting the slepton soft mass which is assumed to be 
(approximately) flavor-independent. There are also the contribution from the finite 1-loop 
graph involving squark or slepton exchange in the /!» = basis: 

A ilm A jnk Y lk* ( H l) {Atnim) + ^(H* 



m 



loop 



16tt 2 



m , 



(16) 



+ 



A iam A< jf3k A ^*ak(^) (Q/3m(-^) + ^Offmi^) 



171- 
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where the Greek indices (a,/3, ...) run from to 3 with Lq = Hi, while the Roman indices 
run from 1 to 3. Here K% k = -A e ojk = Yg, Cf ok = -C* jk = A\ k , and m\ and mf 
denote the squark and slepton soft masses which are assumed to be (approximately) flavor- 
independent. Note that all the parameters in ([14]) and (|16"D are renormalized at the weak 
scale. The contribution involving the sneutrino VEV (Li) = (z/j) in the loop mass has been 
overlooked so far, but turns out to be crucial to fit 7712/^3 ~4x 10 -2 in our framework. 

If U{l)x is anomalous, which is the most interesting possibility, the quadratically di- 
vergent Fayet-Iliopoulos coefficient X 2 Mp naturally yields (4>)/Mp m A. It also leads to a 



nonvanishing U{l)x -D-term [|lq , |l6 |: 



D x « \F\ 2 /M 2 X « \F\ 2 /g 2 \ 2 M^ 



(17) 



where F denotes the SUSY-breaking F-term, Mx ~ g(4>) the U(l)x gauge boson mass. 
The soft scalar mass of $ 7 then receives a D-term contribution 5mj = X(&)Dx- In 
gravity-mediated SUSY-breaking models, this D-term contribution dominates over the stan- 
dard D-term contribution \F\ 2 /Mp. If the U(l)x charge X($ 7 ) is flavor-independent, the 
scalar masses dominated by the D-term contribution would be (approximately) degenerate, 
thereby avoid the dangerous flavor violation [pl|[l7]]. However in our framework, X($ 7 ) are 
flavor- dependent to explain the fermion mass hierarchy. When the D-term contribution is 
important, the requirement to avoid dangerous flavor violation while explaining the quark 
and charged lepton mass spectrum through flavor-dependent severely constrains the 



possible U(l)x charge assignment |18| , and actually leads to the so-called "more" mini- 
mal supersymmetry |L9|| . However the resulting X($ 7 ) do not fit in with our framework 
explaining the small B/L- violating couplings by means of U(l)x- It thus appears that 
gravity-mediated models with U(l)x do not fit in well with our framework. 

The bothersome flavor- dependent U(l)x D-term contribution becomes negligible in 
gauge-mediated SUSY-breaking models with a messenger scale M m <C j^Mx for which 



\JDx "C m so f t « j-F/M m . To discuss the neutrino mass matrix in gauge-mediated case, 
let fi a L a H 2 and B a L a H 2 denote the bilinear terms in the superpotential and soft scalar 
potential in generic basis, and m 2 a ^L a L*^ the soft scalar masses of L a = (Hi, Li). If the 
messenger gauge interactions do not distinguish Hi from Li, B a is naturally aligned to \i a 
and also = m\b a p at M m . In this case, Dj and rn\ iH ^ can be simultaneously rotated 
away as fii at M m , i.e. Dj(M m ) = m 2 LiHi (M m ) = in the basis of fii = 0, and their low 
energy values at Mz are determined by the RG evolution which is governed by the AL = 1 
Yukawa couplings A^ e k and generic L-conserving couplings. The L-violating trilinear soft 
scalar couplings {Cf-D at M z are also determined by the RG evolution. Thus in gauge- 
mediated models, all renormalizable L-violating couplings at Mz are calculable in terms of 
Aijl and also of generic L-conserving couplings. 

Soft parameters in gauge-mediated models typically satisfy: M a /a a ~ mq/a.3 ~ 



S 



rrij/aip at the gauge messenger scale M m where M a , rriq, and m ; ~ denote the gaugino, squark 
and slepton masses, respectively, and a a = g 2 /47r for the (SU(5)-normalized) standard 
model gauge coupling constants g a (a = 1,2,3). Trilinear scalar coefficients do vanish at 
M m and thus their low energy values are determined by the RG evolution. The size of the 
bilinear term BHiH 2 in the scalar potential depends upon how fi is generated. An attractive 
possibility in this regard is B(M m ) = for which all CP- violating phases in soft parameters 



at Mz are automatically small enough to avoid a too large electric dipole moment |22| , |23 
In this case, the RG-induced low energy value of B yields a large tan (3 ps (m% 1 + m 2 H2 + 
2/i 2 )/ BiMz) = 40 ~ 60 which corresponds to x = in view of tan/5 ~ X x mt/mb- In fact, 
a careful analysis of the neutrino mass matrix implies that when x > 1 it is rather difficult 
to fit 1712/m^ «4 x 10~ 2 for reasonable range of soft parameters in gauge-mediated models 
\% without a sizable cancellation [^TJ, and thus here we concentrate on x = 0. 



Analyzing the neutrino masses (|T4|) and (|T^) determined by the RG evolution of couplings 
with the boundary conditions that trilinear soft scalar couplings, B, Bi and , m 2 L i H 1 are a U 
vanishing at M m , and also M a /a a ~ m^/o^ ps mj/ai^ at M m , it is straightforward to find 
that (for x = and thus tan f3 = 40 ~ 60) 

mf* = lO-^tVa,- {^^j , (18) 

where Oj ~ ^A^ 33 w A Zi_?l1 , t = ln(M m /m;-)/ln(10 3 ) and ^1 is the coefficient of order unity 
summarizing the uncertainty of our estimate. Among various terms in the loop mass (|16|). 
the leading contribution to the loop mass comes from the piece involving (Lj)(if 2 ) for large 
tan/?, because Y e (Li)/A e (Hi) ps tan j3 > I. We then have 

m£* « 10- 2 6t 2 nn 3 A^(^A -33 + (^) , (19) 

where the smaller contributions are ignored and again the coefficient £2 of order unity is 
introduced to take into account the uncertainty of our estimate. Here the powers of t oc 
\n(M m /mj) are from (pi) oc Bi(mf) oc t 2 under the boundary condition Bi(M m ) = 0. The 
above neutrino mass matrices are derived in the basis for which Yf A and are diagonal, 
and Yb = Y% 3 and Y T = Y% 3 . At any rate, from (18) and (p!9|), we find 
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rn 



tree 
23 




where we have used Y,- ~ tan/3/95, ~ tan/3/50. Since A233 and A233 are comparable to 
each other, the above result shows that mfj ee gives a dominant contribution. 

Obviously m\j ee is a rank 1 matrix, and thus the total neutrino mass matrix takes the 
form (^) with an approximately singular matrix when m*J ee dominates. We then find 
from (|18D and the following mass hierarchies: 

m 3 « U&rig»U jZ « 10- 1 r / M z A 2fe - hl) , 

loop 

_ 7-T looprr _ m 23 

'''233 

m a « Uam^Uj! « m 2 A 4 , (21) 

where 77 = £i(ln ^/ In 10 3 ) 4 (M z /i 2 /m|). For m z ~ « 200 ~ 400 GeV and /i « 2m,- which has 
been suggested to be the best parameter range for correct electroweak symmetry breaking 
23fl, 7] is roughly of order unity and then the experimentally favored m 3 ^5x 10 -2 eV can 



be obtained for 

7</ 3 -/>i<9. (22) 

Note that in our framework small m 2 /m3 is essentially due to the loop to tree mass ratio, 
while the other small mass ratios mi/m 2 ~ A 4 and m^/Mz ~ 10 — 1 A 2 ^ 3- are from the 
U(l)x selection rule. 

We found many possible U(l)x charge assignments producing all fermion masses and 



mixing discussed so far, while satisfying all the bounds on B/L- violating couplings [24 
through the U(l)x selection rule under the condition that the maximum U(l)x charge is 
not unreasonably large for X(X) = —1. For more detailed discussions, see pif . In this 



paper, we pick two representative solutions: Model 1 and Model 2 which are listed in Table 
1. 

To conclude, we have studied the neutrino mass matrix in supersymmetric models in 
which the observed quark and charged lepton masses and also the suppression of B/L 
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violating couplings are all explained by horizontal U(l)x symmetry. A particular attention 
was paid for the possibility that the neutrino masses and mixing angles suggested by recent 
atmospheric and solar neutrino experiments arise naturally in this framework. It is found 
that our framework fits in best with gauge-mediated SUSY breaking models with large 
tan P ps 50, and favors the small angle MSW oscillation of solar neutrinos over the large angle 
just-so oscillation. Combining the informations from neutrino oscillation experiments with 
those from the quark and charged lepton sector and also the constraints on B/L- violating 
couplings, we find the U(l)x charge assignments producing all the fermion masses and 
mixing angles correctly. This framework determines the order of magnitudes of the neutrino 
mixing matrix elements and mass eigenvalues to be: [/ ps £/ T with U e2 ps U e3 ps A^ 3 Zl = A 2 , 
ps A' 3 "' 2 ps 1 and m l /m 2 ps A 2 ^ 3- ' 1 ) ps A 4 , m 2 /m 3 ps (Loop/Tree) ps £>(10~ 2 ) for 
m 3 ps Ofdir 1 ) x M z \ 2 ^- hl ^ ps 5 x 10~ 2 eV. 
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work is supported in part by KOSEF Grant 981-0201-004-2, MOE Basic Science Institute 
Program BSRI-98-2434, KOSEF through CTP of Seoul National University, and KRF under 
the Distinguished Scholar Exchange Program. 
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TABLES 



TABLE I. U(l)x charge assignments for the MSSM fields of two representative models. 



Model 




Ui,U2,U 3 


di,d 2 ,d 3 


hi hi h 


ei,e 2 ,e 3 


hi,h 2 


1 


8,7,5 


-3,-6,-8 


-1,-2,-2 


7,5,5 


1,0,-2 


-3,3 


2 


7/2,5/2,1/2 


-1/2,-7/2,-11/2 


11/2,9/2,9/2 


5,3,3 


5,4,2 


-5,5 
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